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Abstract. 

Background Neutrino-induced pion production can give important informationon the axial coupling to nucleon resonances. 
Furthermore, pion production represents a major background to quasielastic-like events, one pion production data from 
the MiniBooNE in charged current neutrino scattering in mineral oil appeared higher than expected within conventional 
theoretical approaches. 

Purpose We aim to investigate which model parameters affect the calculated cross section and how they do this. 

Method The Giessen Boltzmann-Uehling-Uhlenbeck (GiBUU) model is used for an investigation of neutrino-nucleus 
reactions. 

Results Presented are integrated and differential cross sections for l;r+ and l;r production before and after final state 
interactions in comparison with the MiniBooNE data. 

Conclusions For the MiniBooNE flux all processes (QE, l;r-background, A, higher resonance production, DIS) contribute to 
the observed final state with one pion of a given charge. The uncertainty in elementary pion production cross sections 
leads to a corresponding uncertainty in the nuclear cross sections. Final state interactions change the shape of the muon- 
related observables only slightly, but they significantly change the shape of pion distributions. 

Keywords: neutrino reactions, nuclear effects, pion production 
PACS: 13.15.-i-g, 25.30.Pt 



INTRODUCTION 

At Nuint 2009 a broad comparison of generator predic- 
tions for QE and one pion cross sections [1] has shown 
major discrepancies (« 30% for QE, up to 100% for pion 
production) between the predictions of various genera- 
tors presently being used by neutrino experiments. At 
this conference (see the contribution of P. Rodrigues), 
a similar comparison shows more consistent results, be- 
cause now each neutrino event generator can be fine 
tuned to the recent MiniBooNE data on charged [2] and 
neutral [3] pion production in charged current (CC) neu- 
trino scattering. 

Does the fine tuning reveal the physics behind it and 
contribute to our understanding of the microscopic mech- 
anism inside nucleus? Is this fine tuning done consis- 
tently for pion production and QE scattering? Answering 
these questions requires that not only the results (as have 
been done at this conference), but also the physical mod- 
els behind the generators are compared. One would need 
to identify the crucial physical inputs responsible for the 
key features of predicted cross sections and to use those 
inputs consistently. We hope that a fruitful discussion be- 
tween different generator-building groups is coming and 
this discussion starts here. 



INPUT 

Accuracy of elementary input determines 
the accuracy of calculated cross sections 

A necessary input into all calculations of neutrino- 
nucleus cross sections are the elementary reaction cross 
sections. For the one pion production the most relevant 
ones are the ANL [6, 5] and BNL [4] data for scatter- 
ing off hydrogen and deuterium. Until our recent publi- 
cation [7] we had used the ANL data as default input in 
the GiBUU code, because for these the absolute values of 
the cross section for da/dg^ were given. In this talk we 
now also show results obtained by using the BNL pion 
production cross sections as input in order to explore the 
consistency of both input data sets with the recent Mini- 
BooNE data and to explore the sensitivity of the nuclear 
data to the elementary interactions. 

Fig. 1 shows that the data obtained in the BNL ex- 
periment are generally higher than those from the ANL 
experiment. In order to obtain an estimate for the sys- 
tematic uncertainty in our calculations, we fit here both 
of these data sets separately. The curves in Fig. 1 show 
the theoretical description, obtained from these fits and 
used as input in GiBUU. We consider them as lower and 
upper bounds of the elementary cross section. 






FIGURE 1. Single-pion production cross section on proton and neutron targets obtained in the BNL [4] (circles; solid lines) and 
the ANL experiments [5] (diamonds; dashed lines). The curves give the lower (ANL-tuned) and upper (BNL-tuned) boundaries on 
the elementary input as used in GiBUU. 



How medium modifications of the A 

resonance properties influence the 

calculated cross section 

The GiBUU simulation code, which besides describ- 
ing many other nuclear reactions can also be used as neu- 
trino event generator was already described at the talk at 
this conference. There we have also emphasized the im- 
portance of final state interactions of outgoing hadrons 
during their propagation throughout the nucleus. The 
verification was provided by comparison with the photo- 
production data for neutral pions. Here we use the same 
set of data to demonstrate the importance of taking into 
account the medium modification of A resonance proper- 
ties. 

In nuclear medium the width F of any baryon acquires 
a collisional contribution 
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FIGURE 2. Momentum distributions of the outgoing pions 
for inclusive K production in scattering of photons of energies 
0.2 -0.8 GeV carbon. 



where the selfenergy Ecoii is a function of density and 
momentum. For the A resonance the model for Ecoii 
was proposed by Oset and Salcedo (OS) in [8]. All 
calculations in this paper, except otherwise noted, have 
been done using this collisional width. 

Fig. 2 compares the photoproduction data with the re- 
sults of GiBUU calculations obtained with and without 
OS medium modifications. Calculations using the OS 
model (red solid lines) are in general agreement with the 
data [9] (red circles), while calculations without medium 
modifications (blue short-dashed lines) noticeably over- 
estimate the cross section in the A peak region. 

The influence of A medium modifications on neutrino 
reactions is illustrated in Figs. 3, 4. The elementary input 
is the cross section for 6 free neutrons and 6 free protons, 
all at rest (cyan dash-dotted line labeled "6p-i-6n"). The 
nuclear effects, namely the Fermi motion, the binding 
nuclear potential and the Pauh blocking of the outgoing 
proton, decrease the cross section by around 5% (red 
solid fines labeled "'^C"). 



Taking into account the OS modification of Delta 
properties (blued dashed curve labeled "'^C, OS") leads 
to an additional decrease of the A-production cross sec- 
tion by another 5 — 8%. This decrease is a consequence 
of the shift of strength from the A peak position to 
larger masses due to the increased broadening. This ex- 
tra strength at the higher masses is, however, cutoff by 
the form factors so that the net effect is a lowering of the 
cross section. 

A similar decrease, as expected, is observed for the 
one pion cross section, with the effect being even larger: 
the cross section is decreased by 15 — 20%. This ad- 
ditional suppression is due to the fact that the A now 
has an increased collision width so that it can undergo 
A + A^ — > A^A^ before the pion decay has taken place. The 
same process, on the other hand, increases the cross sec- 
tion for events with no pions in the final state. We stress 
that it is essential to treat the collisional broadening of the 
A spectral function consistently with the AN collisions. 

Fig. 4 shows the kinetic energy energy distribution 
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FIGURE 3. Cross sections for A production, 1 ;r-production 
originating from A and O;r-production originating from A with 
and without the OS medium modification [8]. 
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FIGURE 4. Kinetic energy distribution of n^ produced in 
neutrino scattering off carbon through the weak production of 
the A resonance and its following decay. The neutrino energy 
is iiv = 1 GeV. The curves labeled OS were obtained using the 
OS in-medium coUisional width of the A. 



of the 7t^ originating from A decay. Before FSI the 
distribution has a rather broad peak at T w 0.1 GeV, 
which is decreased by about 5%, if the OS modification 
of the A properties is taken into account. 

After FSI, the OS modification of the A properties 
decreases the pion distribution significantly more, by 
more than 20% at the peak. This is due to the fact 
that together with the collisional broadening of the A 
resonance explicit 2-body and 3-body collision terms are 
now active. Thus, before having a chance to decay to a 
pion a significant part of As is absorbed in the nucleus. 
The overall effect is a decreased pion production cross 
section which increases the number of pionless events 
and thus serves as a source of the fake QE-like events. 
It is worthwhile to point out already here that the special 
shape in the pion momentum distribution after FSI is not 
affected by the OS modification. 



OUTPUT 

In the following discussions the data from [2, 3] are plot- 
ted versus reconstructed neutrino energy, while the theo- 
retical curves presented here are versus real energies. We 
have found that the energy reconstruction method used 
in the experiment [3], that relies only on the kinemat- 
ics of the outgoing muon and pion, is quite reliable, so 
we neglect the difference. For the later comparisons it 
is also essential to note that all the experimental cross 
sections for ;r+ production were obtained with the full 
MiniBooNE flux. Thus, for positively charged pions the 
distributions with respect to muon or pion energy do not 
depend on any energy reconstruction scheme. This is not 
so for the n^ data where a cut for the (reconstructed) neu- 
trino energies was imposed: only neutrino energies be- 
tween 0.5 and 2.0 GeV were taken into account [3]. Since 
this cut is not possible without a generator-based recon- 
struction procedure the tt" data may thus contain some 
model dependence. We also note that all our earlier cal- 
culations published in Refs. [10, 1 1, 12, 13] used the full 
MiniBooNE flux, without this cut, for all charge states. 
The net result of using now the neutrino energy-window 
for tt" is an increase of the calculated flux-averaged cross 
sections for neutral pions compared to the earlier calcu- 
lations. 



Inclusive cross section and lepton 
observables 

The cross sections for one pion production versus 
neutrino energy for CC l7r+ and l;r" production are 
shown in Fig. 5. 

Figs. 5a,b,c show the results for l7r+ production. As 
expected, for a nuclear target (CH2 in this case) the BNL 
elementary input leads to a significantly larger cross sec- 
tion than the ANL input, simply reflecting that fact that 
already the elementary BNL cross sections are about 
30% larger than those from the ANL experiment (see 
Fig. 1). For both inputs the calculated cross section after 
FSI (red solid lines) is about 30% lower than the cor- 
responding cross section before FSI (blue short-dashed 
lines). 

Before FSI, the most of 1 n^ events come from initial 
A resonance production and its following decay (cyan 
dash-dotted line). This channel is dominant up to about 
0.9 GeV incoming neutiino energy. Some events are 
background ones (dotted line in Fig. 5b). There is also a 
small contribution from higher resonances, and at Ey > 
1 .2 GeV DIS processes start to play an increasing role. 

FSI noticeably decrease the A-originated one pion pro- 
duction (as comparison of panels (b) and (c) shows) due 
to the absorption A^A -^ NN; the similar process is pos- 
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FIGURE 5. (Color online) Integrated cross section for l;r+ (a,b,c) and l;r (d,e,f) CC production versus neutrino energy. Panels 
(a) and (d) contain the results of calculations using ANL- (lower curves) and the BNL-tuned (upper curves) elementary input both 
before and after FSI. Panels (b) and (e) show the contributions of various event origins to the calculated cross sections before FSI. 
Panels (b) and (f) show the same after FSI. Here the BNL-tuned input has been used. Data are from [2, 3]. 



sible also for other resonances. Once a pion is produced, 
independent of its origin, it may also undergo a charge- 
exchange Ti^n -^ Ti^p process, which depletes the ;r+ 
channel as the dominant one. The ensuing reduction of 
the In: -background channel due to FSI is also notice- 
able. Other possibilities for pions to disappear include 
nN -^ coN, ^N, T.K, AK. 

A minor amount of pions comes from the initial QE 
vertex (green long-dashed line), which is only possible 
due to FSI, when the outgoing proton is rescattered. Here 
the main contribution is from the pN -^ N'A -^ N'N n 
reaction. Other possibilities to create pion during the FSI 
would be (oN — > nN, (j)N — > TtN, nN — > nnN. 

For the Itt" production in Figs. 5 d,e,f the curves 
before and after FSI are not so different. This is mainly 
due to the fact, that side feeding from n^n — > n^p, 
which decreases the charged pion output, simultaneously 
increases the neutral pion output. The reverse process 
gives only minor relative contribution, because the initial 
n'^ production cross section is around 5 times lower. This 
is, however, partly compensated by charge exchange to 



the n^ channel through n^n -^ n^ p and other channels 
mentioned above for charged pion production. 

For the 1 ;r'', the A channel is less important than in the 
case of charged pions, while background processes, and 
at higher energies higher resonances and DIS play a rela- 
tively larger role than for 7r+. These two latter processes 
ensure the change in slope at about 1 .2 GeV, which re- 
duces the the discrepancy between the theoretical curves 
and data. 

The MiniBooNE collaboration has recently also pub- 
lished distributions versus muon kinetic energy, da/dT^, 
muon angle with respect to the neutrino beam direc- 
tion, da/dcosOf^, and squared four-momentum trans- 
fer, da/dQ^ [2]. All of them are calculated within the 
GiBUU and compared to the data in Ref. [7]. Here we 
present only the kinetic energy distributions in Fig. 6. 

It shows that FSI hardly influence the shape of the 
distributions versus muon kinetic energy, the same is 
true for other muon observables. This insensitivity is due 
to the fact that the only effect of FSI on these muon 
observables is that they can remove events in which an 
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(Color online) Kinetic energy distributions of the outgoing muons for the 1 ;r+ and 1 ;r production in the MiniBooNE. 
[2, 3]. The blue dashed curves give the results before FSI, the red solid curves those with all FSI included. 



initially produced pion (or A) was later on reabsorbed 
and bring in events in which the pion was produced only 
during FSI. 

For 1 K^ production, when referring to the shape-only 
comparison, the calculated curves correspond to the data. 
However, the muon kinetic energy distribution is consis- 
tently 20% (upper boundary) to 60% (lower boundary) 
lower than the data, see Fig. 6. For Itt" production, the 
shape of the curves evidently differs from the shape of 
the data. 



Pion Observables 

Fig. 7 presents the results of our calculations for the 
kinetic energy distribution of the outgoing pions. Unlike 
the case of muon-related observables, for pion-related 
distributions the FSI noticeably change their shape. Sur- 
prisingly, the shape before FSI is similar to that observed 
in the data. Nevertheless, the final state interactions must 
be inevitable for nuclear targets. After FSI, however, the 
shape of the calculated distributions looks markedly dif- 
ferent. In particular, there is a significant lowering around 
T^ w 0.2 GeV for ;r+ and at around pi w 0.3 0.3 GeV, as 
a direct consequence of the A;rA^ dynamics in nuclei. 

The following processes are important for this struc- 
ture: pion elastic scattering in the FSI decreases the pion 
energy, thus depleting spectra at higher energies and ac- 
cumulating strength at lower energies. Simultaneously, 
there is charge exchange scattering. At the same time 
pions are mainly absorbed via the A resonance, that is 
through TtN — ^ A followed by AA^ — > A^A^, which leads 
to the reduction in the region of pion kinetic energy 
0.1 —0.3 GeV. For tt" production the additional increase 
of the cross section at lower energies comes from the side 
feeding of the n^ channel from the dominant 7r+ channel 
due to the charge exchange scattering n^n 



ny. 



Inverse feeding %^p -^ n^n is suppressed, because at 
the energies under consideration, about 5 times less ;r''s 
than ;r+s are produced. The change of the shape of the 
spectra due to FSI is similar to that calculated for neutral 
current Itt" production in [10, 11]. 

As we already emphasized in another talk at this 
conference, the particular shape calculated here for the 
neutrino-induced pions is in line with that observed ex- 
perimentally in (7, tt") production on nuclear targets [9]. 
It has also been observed in the independent calculations 
of Hernandez et al. [14]. Since the shape depends on FSI 
and since the FSI are the same in both neutrino-induced 
and photon-induced reactions the absence of this special 
shape in the neutrino data is surprising. 

Fig. 8 shows the pion angular distribution of the of 
the CC 1 n^ and 1 n^ production. Here data are available 
for tt" only; for forward scattering they are significantly 
higher than the lower boundary of our calculations, and 
only slightly higher compared to the upper boundary. 
The forward peaking is reproduced, but it is not as strong 
as exhibited by the data. 



CONCLUSIONS 

In conclusion we summarize here a few essential points. 

• Any neutrino event generator should have all ini- 
tial reaction mechanisms (quasi-elastic scattering, A 
production, production of higher resonances, back- 
ground 1- and many -pion production, DIS) under 
control. Even the MiniBooNE data require the in- 
clusion of processes beyond the A excitation. 

• There is a strong dependence of theoretical results 
on in-medium Delta properties. 

Finally, we stress that because of the strong depen- 
dence of the results on the elementary pion production 
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FIGURE 7. (Color online) Kinetic energy distribution of the outgoing ;r+ and momentum distribution of the outgoing ;r for one 
pion production at MiniBooNE. Data are from [2, 3]. The curves are as in Fig. 6. 
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FIGURE 8. (Color online) Distribution of the outgoing n^ and ;r in their angle relative to neutrino beam for one pion production 
at MiniBooNE. Data are from [3]. The curves are as in Fig. 6. 



cross section new, improved data for neutrino-induced 
pion production on elementary targets are clearly needed 
to resolve the ambiguity in pion production on nuclear 
targets. Only when such data have become available the 
question of a better determination of the axial coupling 
to nucleon resonances can be approached again. 
This work is supported by DFG and BMBF. 
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